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ABSTRACT 

Three point mass gravitational lens equation is a two-dimensional vector equa- 
tion that can be embedded in a tenth order analytic polynomial equation of one 
complex variable, and we can solve the one variable equation on the source tra- 
jectories using recipies for Fortran or C (portable for C++ or Cjj) in Numerical 
Recipes, or using packages such as Mathemetica, Matlab, etc. This ready solv- 
ability renders fitting microlensing light curves including triple lenses a normal 
process, and such was done in a circumbinary planet fit for MACHO-97-BLG- 
41. Subsequently, there was a claim that converting the triple lens equation into 
the analytic equation was rather cumbersome, and the impressionable judgement 
has caused an effect of mysterious impedance around the perfectly tractable lens 
equation. There are judgements. Then, there is nature. We looked up for one 
of the quantities of highest precision measurements: electron g-factor correction 
a e = g/2 — 1. The current best experimental values of a e agree to eight signif- 
icant digits with the theoretical value, and the theoretical calculation involves 
more than one thousand Feynman diagrams - many orders of magnitude messier 
than the triple lens equation coefficients. We seem to have only choice to be 
compliant to nature and its appetite for elegant mess and precision numerics. In 
fact, the triple lens equation coefficients take up less than a page to write out 
and are presented here for users' convenience. 

Subject headings: gravitational lensing: 



1. Gravitational Triple Lens Equation 

The lens equation of a set of three gravitationally bound point masses is written with 
three (real) relative mass parameters, E\, and £3, and three complex parameters for their 
2-dimensional positions. If z is the position of an image and u is the position of its source, 
the lens equation reads as follows. 

u = z — - — X — — — ^— = Z — f(z]Xj) (1) 

Z — X\ Z — X 2 Z — £3 
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The mass parameters are subject to a constraint e\ + e 2 + e 3 = 1 where 1 = M is the 
total mass, and the position variables include coordinate degrees of freedom: two degrees 
of freedom for translation and one degree of freedom for roation in the two-dimensional lens 
plane. We choose a coordinate system such that the position of a mass, xi, and the center 
of mass of the other two elements, X4, define the lens axis along the real axis of the complex 
plane. 

(e 2 + e 3 ) x 4 = e 2 x 2 + e 3 x 3 (2) 

Then, x\ and X4 are real, and the triple lens system is completely specified by a set of five 
parameters. 

£=\xi-x 4 \; £' = \x 2 -x 3 \; 7 = [0, tt] ; e 2 = (0, 1) ; e 3 = (0, 1) (3) 

The angle 7 is the angle between (x 2 — x 3 ) and the lens axis. In order to find the images of 
a given source, we need to solve the lens equation. 

When the lensing system and a source star are far apart, there are four images of the 
source star, which we can verify easily for a source at u = 00. From the lens equation (1), 
we find that the images of the source at uj = 00 are at the position of the source, z = 00, 
and at the three lens positions, z = x±,x 2 , and x 3 . The magnifications of the images at the 
lens positions are 0, and they are hardly images in any practical sense because there are 
no photon fluxes related to the image positions. In other words, there is only one physical 
image for a source at 00 - namely, the image of the unmagnified source star. However, if we 
move the source toward the lens system, the images at the lens positions move away from 
the lens positions, and they do have non-vanishing photon fluxes. Thus, we freely speak of 
image positions with zero magnification as the continuity limit so that we don't have to cut 
out the lens positions from the image space. In practice, 00 is a large distance limit, which 
may be no more than a few hundred au in any physical relevance of Galactic lensing. 

If we consider a source trajectory away from the caustic regions, the images form four 
smooth curves anchored at the base points of z — x±,x 2 ,x 3 , and 00 which are the image 
positions of the source at u> = 00. The complex lens plane can be considered a large two 
sphere with one point at 00, and then, the four smooth disjoint image curves are loops with 
fixed points at the base points. The image curves can be calculated esaily from the base points 
using Newtonian transportation method (the algorithm can be found in Numerical Recipes). 
However, the physical interest of lensing lies in caustic regions where the lensing signals 
are most obvious, and this makes the Newtonian method useless in microlensing business. 
The large lensing signals arise where the inverse Jacobian determinant of the lens equation is 
large, and the large Jacobian determinant makes Newtonian transportation method unstable 
and unusable. The method of infinitesimal corrections based on derivative values is valid 
where the Jacobian determinant is practically finite and the mapping is non-degenerate. 
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Laguerre's method (see Numerical Recipes) is based on complex analysis and uses first and 
second derivatives of the polynomial to help converge to the solutions from initial values 
which one is free to guess. 

Another complication of the caustic region is that the size (which may be 1 — 10/ias) 
of the source star can not be ignored when the source star crosses a caustic curve. In fact, 
the luminosity profile of the source star can manifestly affect the shape of the light curve, 
and the luminosity profile dependence of the light curve shape during a line caustic crossing 
lasts ~ 3 stellar radius crossing times (Rhie and Bennett 1999). In order to incorporate the 
finite size effects, one can locally pixelize the image plane (z-plane) in the neighborhood of 
the critical curve near the images and count up the image pixels that are mapped into the 
source disk (Bennett and Rhie 1996). The luminosity profile of the source can be easily 
incorporated by weighting each image pixel by the luminosity shape function value at the 
corresponding source position. The pixel size can be adjusted for a desired resolution. In the 
ray shooting method which has been designed mainly to handle a system of a large number 
of lensing elements solving whose lens equation is impractical if not impossible, the entire 
lens plane is pre-pixelized and each pixel in the image plane is tested whether the center of 
the pixel is mapped into the area defined by the source disk in the source plane. The ray 
shooting method is inefficient for fitting light curves of low-multiplicity point mass lenses as 
microlensing planet systems. 

The astro-ph version of Gaudi, Naber, and Sackett (1998) conveys an impression that 
solving the lens equation by finding the roots of the 10-th order polynomial analytic equation 
is slower than using the ray shooting method for triple lens light curve fitting. However, 
Gaudi (private communication) recently informed us that the comment on the calculation 
speed was a comparison between binary lens and triple lens but not a comparison between 
root finding method and ray shooting method for the triple lens equation. 

2. The Tenth Order Polynomial Equation 

The lens equation is an explicit function from an image position z to its source position 
u, and u(z,z) is a genuine real function, namely a function of both z and z. However, the 
z-depdence is linear, and this simplicity is behind the analyticity of the differential behavior 
which can be completely described by one analytic function k = d z uj. The linearity in z 
also makes it easy to find an analytic equation where the lens equation is embedded. Using 
z — uj = f(z; Xj) and z — cu = f(z; Xj), 



z-u = f(f(z;xj)+u;xj) 



(4) 
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If we let H = ziz 2 z 3 and G = ei^^ + 62^3-21 + £3-21-22 = J2 cyc e i z j z k, then f = G/H, and it is 
simple to see that equation (4) is a tenth order polynomial equation. If we let cuj = Cu — Xj, 



= (z-cu)(G + iDiH)(G + u 2 H)(G + u 3 H) - ^ He t (G + UjH)(G + u k H) (5) 



An analytic polynomial equation has the same number of solutions as the order (see any 
textbook on complex variable or mathematical physics), there can be up to ten images in a 
triple lensing. There are only four images for uj = 00, the number of images changes by two 
at a caustic crossing, and the caustic curves form heierachical structures of domains for high 
multiplicity images. There are triple lenses with domain T> 3 the sources therein produce ten 
images, the maximum possible number of images (Rhie 1997). As we repeatedly emphasized 
in Rhie (1997), the triple lens equation is equivalent to the tenth order polynomial analytic 
equation only in the domain V 3 , and the statement on the equivalence in section 3 of Gaudi, 
Naber, and Sackett (1998) seems to be a misquote which can mislead the readers to think 
that the two equations are equivalent everywhere. 

In order to fit a triple microlensing light curve, we need to solve the tenth order analytic 
equation (and select the image solutions that satisfy the lens equation). This can be done 
numerically using root finders available in the literature, and we only need to type in the 
coefficients. The coefficients may appear to be cumbersome as declared in Gaudi, Naber, 
and Sackett (1998), and it is indeed the case if we, for example, calculate the coeffients using 
Mathematica. The output of an algebraic computing package is (unnecessarily) messy even 
for the binary lens equation. Thus, it is useful to group (or not to unfold) the coefficients, 
which is a natural intermediate process in hand calculations. 

In the center of mass system, ejXj = 0, H and G are written with four coefficient 
functions, a, b, c, and d: H = z 3 + az 2 + bz + c and G = z 2 + az + d, where a = —{x\ + x 2 + 
£ 3 ), b = xix 2 + X1X3 +x 2 x 3 , c = -xix 2 x 3) and d = Y^ cyc e i x j x k- If we let a w = a>i + lj 2 + ^3, 



b u = UJ1U2 + + c w = U1U2UJ3, and d w = J2 cyc eiUjUk, then equation (5) becomes 

= G\z -u) + G 2 H((z - c^K + 1) + GH 2 ((z - u)b w + au) -u)+ H 3 ((z - u) Cu) + dj) (6) 

If we let G 3 = H ok z k (k < 6), G 2 H = H lk z k (k < 7), GH 2 = H 2k z k (k < 8), and 
H 3 = H 3k z k (k < 9), the equation becomes 



eye 



10 




(7) 



k=i 



where the polynimial coefficients are 



• eff (k) = ( f/ofc-i + i^ife-i a w + H 2k ~i b u + H 3k ^ cj) 

- (H ok u + H lk (ua^ - 1) + H 2k (ub^ + a w - u) + H 3k (wc w + b u )) 
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The coefficients Hy are polynomials of a, b, c, and d where the polynomial coefficients are 
simple combinatoric integers. 

• # 39 = 1; # 38 = 3a; # 37 = 3b + 3a 2 ; # 36 = 3c+Qab + a 3 ; H 35 = Qac+3b 2 + 3a 2 b; # 34 = 
6bc + 3a 2 c + 3ab 2 ; H 33 = 3c 2 + 6abc + 6 3 ; H 32 = 3ac 2 + 3b 2 c; H 31 = 36c 2 ; H 30 = c 3 . 

• # 28 = 1; #27 = 3a; # 26 = d + 2b + 3a 2 ; # 25 = lad + 4ab + a 3 + 2c; # 24 = 2db + da 2 + 
4ac+2a 2 b + b 2 ; # 23 = 2dc + 2dab + 2a 2 c + ab 2 + 2bc; # 22 = 2cad + db 2 + 2abc + c 2 ; # 2i = 
2bcd + ac 2 ; # 2 o = c 2 d 

• #17 = 1; #i6 = 3a; # i5 = 2d + 3a 2 + b; # M = 4ad + a 3 + 2ab + c; # i3 = d 2 + 2a 2 d + 
2bd + ba 2 + 2ac; # i2 = ad 2 + 2abd + 2cd + ca 2 ; H u = bd 2 + 2acd; H 10 = cd 2 

• #06 = 1; #05 = 3a; # 04 = 3d + 3a 2 ; # 03 = 6ad + a 3 ; # 02 = 3d 2 + 3a 2 d 2 ; H 01 = 
3ad 2 ; # 00 = d 3 



2.1. Comments 

An interested party may download the source file of this manuscript to avoid tying the 
coefficients. We encourage to check the coefficients, however. It is a quick exercise once 
one adopts the poor person's calculation with pencil and paper as shown above; also with 
the free biocomputer which is harder to hack either internally or externally barring the long 
term process of brainwashing. We also found it useful to test the symmetric cases in rh97 
whose image solutions behavior (for example, the number of images) is known. The critical 
curve is obtained by solving k = e l2ip which is an eighth order polynomial equation. The 
caustic curve is obtained by applying the lens equation to the critical curve solution. It is 
fine to use tp — [0, 2tt) as the parameter for equal interval sampling. Let 5 = 2n/N for iV not 
too large and observe the intervals (or speeds) of the solutions on the critical curve and the 
caustic curve. Note especially the density of the solutions around the cusps of the caustic 
curve. It is worth pausing for a moment counting the relative numbers of the solutions on 
the stellar caustic and planetary caustics for planet systems lenses. 

For the measurements and theoretical calculations of the anomalous magnetic moment 
of the electron, we have consulted Peskin and Schroeder (1999) and Kinoshita (1990). We 
have considered drawing all the Feynman diagrams to lay out the degree of lengthy squiggly 
messiness but given up, and our misadventure may be considered an indirect testimony of 
the degree of mud wrestling the exquisite anomalous magnetic moment requires. 

Barring the notion that lensing community may have been chosen to be dealt with laxed 
scrutiny of nature, we have no doubt that we have no luxury to complain about a bit of 
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algebra we encounter in lensing. In fact, we find it a cherished treasure that low multiplicity 
point lenses are exactly solvable and their light curves can be reconstructed and interpreted 
without ambiguities. Microlensing events do share the transiency with the scattering events 
in accelerator particle physics. The both need high resolution data for minute rare prized 
signals and complete data to interpret them. The both rely on methodical analyses building 
from the simpler and dominant events to more rare events. Thus, it is important to have 
a homogeneous and comprehensive data set of microlensing events where consistencies can 
be tested within in order to find microlensing planets. The best bet for such data set is 
microlensing from space. Accelerator particle physics of the last century is at the foundation 
of the Standard Model or the Theory of Matter. We expect that space microlensing planet 
search will lay a foundation of extrasolar planet physics within a few years of operation of a 
small space telescope (Bennett and Rhie 2000). In the comprehensive data set, the so-called 
high magnification events with only stellar caustic signals will form an independent data set 
that can be used for a consistency check of the interpretation of the planetary light curves 
as a subset. 

This note is based on the work with D. Bennett for Bennett et al. (1999). 
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